Noisy Control, the Adiabatic Geometric Phase, and Destruction of the Efficiency of 

Geometric Quantum Computation 
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We examine the adiabatic dynamics of a quantum system coupled to a noisy classical control 
field. A stochastic phase shift is shown to arise in the off-diagonal elements of the system's density 
matrix which can cause decoherence. We derive the condition for onset of decoherence, and identify 
the noise properties that drive decoherence. We show how this decoherence mechanism causes: (1) 
a dephasing of the observable consequences of the adiabatic geometric phase; and (2) the loss of 
computational efficiency of the Shor algorithm when run on a sufficiently noisy geometric quantum 
computer. 
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As the adiabaticgeometric phase approaches its twen- 
tieth birthday Q, Q it continues to impact research in 
physics and chemistry. The original scenario considered 
a quantum system Q with a discrete, nondegenerate en- 
ergy spectrum whose dynamics was driven by a classical 
set of control parameters B(f). It was shown that when 
B(i) was cycled adiabatically through a loop in parame- 
ter space, and Q was initially prepared in an eigenstate 
|-E(0)) of the initial Hamiltonian H(0), its quantum state 
would return to the initial state l-E'(O)) at the end of 
the cycle, to within a phase factor exp The phase 
shift <p contained a geometric contribution now known 
as the adiabatic geometric phase (AGP) which had been 
discarded in prior treatments of quantum adiabatic dy- 
namics. Since then, this original scenario has been gen- 
eralized in a number of ways: (1) the requirement of 
nondegeneracy was removed |3| ; (2) the adiabatic restric- 
tion was removed (3) the control field B(i) was al- 
lowed to have its own dynamics @; and (4) Q was al- 
lowed to interact with a quantum environment as well 
as with the control field B(t) 0, 0, 0- In this Letter 
we return to the original scenario, though we allow the 
classical control field to contain a noise component B„ (t) 
along with the (deterministic) adiabatically varying com- 
ponent B Q (t): B(t) = B Q (t) + B n (t). We will show that 
B n (£) causes a stochastic phase shift to appear in the off- 
diagonal elements of Q's density matrix which, for suffi- 
ciently strong noise, causes decoherence. We derive the 
condition that determines the onset of decoherence (see 
eq. JBJ) , and identify the noise properties that drive the 
decoherence. We show that this decoherence mechanism 
causes all physical consequences of the AGP to become 
unobservable. 

The AGP has also had an impact on the new research 
area of quantum computation. In 2000, a proposal was 
put forth for a geometric quantum computer (GQC) |9j 
which uses the AGP to encode conditional phase shifts 
into the quantum state of a GQC. It was argued that the 
AGP gave the GQC a modest degree of fault tolerance to 
uncertainty in the control parameters, though the need 



to examine the robustness of a GQC to decoherence was 
also noted. In the second part of this Letter we carry out 
such a decoherence analysis. We show that the above de- 
coherence mechanism arising from noisy control severely 
impacts the performance of a quantum algorithm run on 
a noisy GQC. Specifically, we show how noisy control 
causes the Shor algorithm to loose its computational ef- 
ficiency so that the time to factor an integer (AT) grows 
exponentially with the size of the integer (logiV). We 
establish the condition for the onset of decoherence, and 
show that the success probability for the algorithm is 
smoothly degraded as the noise becomes stronger, reach- 
ing the point at onset where the computational efficiency 
of the Shor algorithm has been completely destroyed. We 
also show that at onset of decoherence, the entanglement 
fidelity of a 2-qubit entangled state is degraded to that of 
a mixed state upon output from a sufficiently noisy con- 
trolled phase gate of the type used in a GQC. This noise- 
induced loss of entanglement confirms that the GQC is no 
longer behaving as a quantum computer. To the best of 
our knowledge, this Letter provides the first quantitative 
demonstration of how the computational efficiency of a 
landmark quantum algorithm such as Shor's is destroyed 
when it is run on a sufficiently noisy GQC. A detailed 
presentation of all arguments will be reported elsewhere 

1. General Analysis: We consider the situation where 
the dynamics of an N-qubit system Q is driven by an ex- 
ternal control field B(f). The control field is composed of 
a deterministic component B a (i) and a noise component 
B„(t). We will be interested in cases where B n (£) exe- 
cutes a cyclic evolution with period T. We take the noise 
component B„(i) to be a stationary stochastic process 
with zero mean (B n (i) = 0) and noise correlation time r c . 
To simplify the notation we assume the noise is isotropic. 
This assumption is not essential, and the following analy- 
sis can easily be repeated without it. For isotropic noise, 

the noise correlation function B!^(t)Bn(t') = Sija 2 f(r). 
Here: (1) Sij is the Kronecker delta; (2) a 1 is the vari- 
ance of the noise {B l n {t)) 2 ; (3) r = * - t'; and (4) /(r) is 
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the normalized noise fluctuation profile with peak value 
/(0) = f and temporal width of order r c . Because of 
the stochastic nature of B„ (t) , each application k of the 
control field contains a different realization of the noise 
B n (t; k). To account for the dynamical effects of noise, it 
is necessary to introduce an ensemble of identical quan- 
tum systems Q k , all prepared in the same initial state 
10(0)), but whose evolution is subject to a different noise 
realization B n (£; k). 

Quite generally, the system Hamiltonian H(t) can be 
written as the sum of two terms: H a {t) and H s (t). Here 
H a (t) is the Hamiltonian in the absence of noise and de- 
pends solely on B a (i). The stochastic interaction term 
H s (t) depends explicitly on the noise component B„(t) 
and its action varies from one application of the control 
field to another. We introduce the instantaneous energy 
eigenvalues E n (t) and eigenstates \E n {t)) of the noiseless 
Hamiltonian H a {t): H a (t)\E n {t)) = E n (t)\E n (t)) . Note 
that the periodicity of B n (t) implies that H a (t + mT) = 
H a (t). As in the original AGP scenario, the instanta- 
neous energies E n (t) are assumed to be discrete and non- 
degenerate for all t of interest so that \E n (t + mT)) = 
\E n (t)}. Introducing the minimum energy level separa- 
tion A = min \E n (t) — E m (t)\, the noiseless dynamics 
will be adiabatic if h/T <C A. Furthermore, as shown in 
Ref. , the dynamics associated with the noise will also 
be adiabatic if %/t c <C A. Both of these conditions are 
assumed in the remainder of this Letter. 

As noted above, each ensemble element Q k sees 
a different noise realization B n (t;fc) added to B Q (t). 
The dynamics of Q k is driven by H(t; k) — 
H a (t) + H s (t;k) which generates the final state 
\ip(t;k)) = U(t f ,t ;k) 



gives: 



exp 



-{i/h)f' dtH(t;k) 



■0(0)), where U(t f ,t ;k) = 
Thus the final state produced 
by the noisy control field varies over the ensemble. Conse- 
quently, calculation of the expectation value of an observ- 
able requires an ensemble average over the noise, along 
with the usual quantum mechanical averaging. 

We divide the time interval (to, tf) into J subinter- 
vals of duration e = (tf — to) /J by introducing inter- 
mediate times tj = t + je (j = 0,...,J). Eventu- 
ally we let J — > oo. The propagator U(tf,to) factors 
into a product of propagators: U(tf,to) = U(J) ■ ■ ■ U(0), 
where U(j) — U(tj,tj-i). One can show 0, EI that 
in the adiabatic limit the matrix elements Ui m (j) — 
(£i(ij)l^(i)|£m(^--i)) are: 



Uhn{j) 
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{E l (t j )-nj l (t j ) + H l J(t j )} + o( 



(i) 



Here ii(tj) = i(Ei(tj)\Ei(tj)) is the time derivative of the 
AGP for energy level I; r Zm (^) = (Ei(t;j)\E m (tj)) is the 
non-adiabatic coupling for levels I and m; and H l s l (tj) = 
{Eiit^Hsitj^Eiitj)). Choosing |0(i o )) = \E k (t )); re- 
constructing U(tf,to) from the U(j); and letting J — > oo 



\iP{t f )) = exp[-i{T a (k) + T s (k)}]\E k (t )). (2) 



Here r o (fc) = (l/ft) J% dt{E k {t) - h%(t)} and T s (k) = 
(1/h) f t ! dtH* k (t). Eq. © shows that noise introduces 
a stochastic phase shift r s (fc) into the final state \ip(tf)) 
which is sensitive to the initial energy eigenstate \Ek{to)) 
through H^ k (t) . 

To bring out the physical consequences of the stochas- 
tic phase T s (k) we consider an initial superposition state 
IV'(to)) = J2k c k\Ek(t Q )) . From the linearity of quantum 
mechanics: \ip(tf)) = ^2 k CkU(tf 7 t )\E k (t )). Informa- 
tion about the phase coherence of the final superposition 
is carried in the off-diagonal elements of the final density 
operator p(tf) = \^(tf))(ip(tf)\. One can show [ljj] that 
the noise averaged density matrix elements are: 



Pkj(tf) = c fe c*exp [-iT a {k,j)]D(k,j), 



(3) 



where r a (fc, j) — T a (k) — T a (j), and the decoherence fac- 
tor D(k,j) is: 



D(k,j) = exp [-iT s (kJ)], 



(4) 



with r s (fc, j) = T s (k) — r s (j). We see that the decoher- 
ence factor D(k,j) is a direct consequence of the noise- 
induced phase shift T s . Often, H s (t) is linear in the noise 
component B n (i) so that H s (t) = — (7?i/2)B n (t) • O. 
Here 7 is the coupling constant and O is a vector opera- 
tor. In the adiabatic limit, r c <C (tf — to) = T. Thus, if 
we partition the integration interval (to, tf) that appears 
in the definition of T s (k, j) into M = T/t c subintervals of 
duration r c , we render T s (k,j) into a sum of uncorrelated 
random variables T™Qs,j): T s (k,j) = J2m=i^Ti^: j)- 
Because the noise is stationary, the T™(k,j) have identi- 
cal probability distributions. If the T™(k,j) are not only 
uncorrelated, but also statistically independent, it follows 
from the Central Limit Theorem that T s (k,j) will have 
a Gaussian probability distribution. In this case, one can 
show that T s (k,j) = and 



r2 /, ■ \ Vl 2 o- 2 T 



(5) 



where: (1) rj cycles of the control field have been applied 
(r?T = T); (2) I kj = J Q T dtdt'6 kj (t)-6 kj (t')f(t-t'), with 
6 kj (t) = OkkW-Otiit), and6 fefe (t) = (E k (t)\6\E k (t)); 
and (3) f(t — t') is the normalized noise fluctuation pro- 
file introduced earlier. Having the mean and variance 
of the Gaussian probability distribution for T s (k,j), we 
can evaluate the noise average in eq. (J3J. The result is 
D(k,j) — exp[— T 2 s (k, j)/2]. One expects that when the 



phase uncertainty yT^(k,j) ~ 2ir, the superposition in 
\ip(tf)) will have been dephased/decohered by the noise. 
Inserting T^(k,j) ~ (27r) 2 into our expression for D(k,j) 
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gives D(k,j) ~ 3 x 10~ 9 so that the off-diagonal ele- 
ments of p(tf) (see eq. ©) are effectively zero and the 
noise has in fact caused an effective collapse of the wave- 
function \ip(tf)). One can show |n| that the noise vari- 
ance <7 2 , the average noise power absorbed by Q per unit 
volume P/V, and the effective bandwidth Aoj of the ab- 
sorbed noise power are related: P/V = a 2 Auj. Using 
this relation, together with eq. (JSJ, one can re-express 



the condition for onset of decoherence \/T 2 (k,j) ~ 2n as 
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(6) 



We consider a number of applications of this general anal- 
ysis in the remainder of this Letter. 

2. Dephasing the AGP: In our first application we re- 
consider the original AGP scenario Q , allowing for noisy 
control. Here Q is a single qubit, and B a (t) precesses 
about the z-axis with period T at an angle 0q. Let C 
denote the contour traced out by the tip of B a (t) in a 
time T. To observe the AGP, the initial state must be a 
superposition: |^(0)) = c+\E + (Q)) +c_|£L(0)}. Eq. © 
gives the matrix elements of the final density matrix. The 
diagonal elements are real-valued, and thus contain no 
information about the AGP. The off-diagonal elements, 
however, do depend on the AGP through r o (+, — ) = 
r o (+) - r o (— ), where L a (±) are defined below eq. (2). 
In an NMR experiment, observation of the AGP is carried 
out by measuring the transverse magnetization whose ex- 
pectation value depends on p-\ (t/) 0, 0|. As the 

analysis of Section I showed, sufficiently noisy control 
causes p^ (tf) to effectively vanish so that all physi- 
cal consequences of the AGP thus become unobservable. 
The stochastic phase shifts L s (±) generated by noisy con- 
trol dephase the final superposition state, reducing it to 
a mixture of the states \E±(0)) in which all AGP ef- 
fects are absent. In the NMR setting, the precession of 
B a (i) is produced by varying the phase of the rf mag- 
netic field H r f(t) = B r f{t)~k. If we imagine the noise is 
due to rf power fluctuations, H s (t) — ~(jh/2)B n (t)cr x 
which identifies O = xcr x . One can show that 
for this type of noise J_| = 4t c T sin 2 6>o . If we de- 
note the static magnetic field that splits the nuclear 
energy levels by Bo = BqZ, it is well-known 0] that 
sin 2 6 = B 2 /(B 2 f + {B - B rf ) 2 ). Using this result for 
7_i in eq. (El gives: 



in eq. 



ril 2 T 
4tt 2 



B 2 f 



Bh + (Bo - B rf ) 2 



'rf 



fells;) < 7 > 



Eq. (J7J identifies the noise properties that impact de- 
coherence, and allows a quantitative assessment of how 
much noise can be tolerated before all AGP effects are de- 
phased by noise. We believe that this is the first demon- 
stration of how noisy control causes a dephasing of the 



AGP, and the first quantitative analysis of the onset con- 
dition for this dephasing mechanism. 

3. Shor Algorithm on a Noisy GQC: We now consider 
how noisy control impacts the performance of a GQC. 
The universal set U of quantum gates used to construct a 
GQC contains the 1-qubit Hadamard gate H and the set 
of all possible 2-qubit controlled-phase gates B(<j)) (tp E 
[0, 2tt)). The action a{B(<f>) on the 2-qubit computational 
basis states (CBS) \xy) is: B((f))\xy) — exp[ixy4>]\xy) , 
where x, y = 0, I. The AGP is used to encode the condi- 
tional phase shift </>. We focus on the operation of B{4>) 
in the presence of noise as it is the only gate in the set 
U that can introduce entanglement into the dynamics of 
a GQC. The conditional phase shift 4> is implemented 
using a four part pulse sequence P = P0P1P2P3, with 
Pq (P3) applied first (last). Here P = Cirx; Pi = C1T2] 
P2 = Po', and P3 = Pi, where C is the cyclic evolution 
of B a (i) introduced in Section 2; C is the time-reverse 
of C; and 7Ti (7r 2 ) is a 7r-pulse applied to the first (sec- 
ond) qubit. By appropriate choice of C , any phase shift 
(j) can be produced. The 1-qubit CBS are: |0) = | £_(*)) 



and |1) = \E + (t)). Let k — (ii, i 2 ) with i\, i 2 



0, 



1; then \E k (t )) = ^(io)) ® \E t2 (t )). The pulse se- 



quence P maps k ^> ki — * k 2 — » k% ^ k^, 



with kj = 



(*i®Sz=i ' m °d 2, i2®J2i=o ' m °d 2)' anc ^ ® ^ s addition 
modulo 2. For |^(*o)) = \Ek(to)), eq. © gives the final 
state \ip(t f )) with r«(fc) = (1/fi) £f =0 r -( fc i); r s (fc) = 
(lA)Ef= r s (fei); and F a (fc ; ) = dt[E h (t) - 

hj kl (t)}; T s (h) = dtH*< k >(t), with t (l) =t + lT 

and tf(l) = to(l) + T. For the initial superposition state 
|^(*o)) = Ck\Ek{to)) + Cj\Ej(to)), the final density ma- 
trix p k j(tj) is given by eq. ©■ The decoherence factor 
D{k,j) is given by eq. @ with T s (k,j) = £f =0 T l s (k, j), 
and T l s [k,j) = jZ$dt[H*>*{t) - H™>(t)]. Because 
r c -C T in the adiabatic limit, the T l s (k,j) are un- 
corrected so that r 2 (fc, j) = £?=o( r U fc .j)) 2 - If the 
r[(k,j) are also statistically independent, one can show 
that D(k,j) = exp[-r 2 (fc,j)/2] and T 2 (k,j) is 
given by eq. with hj ~ > J2i=o^kj- Here Aj nas 
the same form as Ikj, except that Okj{t) — > O l h .{t) — 
Ofeifci(i) — Ojiii(i)- The condition for onset of decoher- 
ence is again T 2 s (k,j) ~ 47r 2 . As an example, imagine 
we input the Bell state (1/V2) [|00) + |ll)j into B(<£). 
Then c fe = c, = 1/V2; fc = (00); and j = (11). 
If we assume again that the noise is due to rf power 
fluctuations, one can show 10] that the noise-averaged 
entanglement fidelity F = (ilj(to)\p(tf)\i/j(t )) = | + 

r.(fc) - r„(j). 

shows that 



|cos r a (fc, j)-P (fc, j), with r a (fc, j) 
Since £>(fc,j) = exp[-T 2 s (k, jj/2], eq. © 
it goes smoothly to zero with increasing noise variance 
a 2 . Thus F — > 1/2 which is the entanglement fidelity 
for an initial state which is a uniform mixture of the 
states j 00) and |11), indicating that entanglement is de- 
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stroyed once the control field driving B(<f>) becomes suf- 
ficiently noisy. One can show [l0| that for this Bell 
state ^2i =Q I l k j — 32T c Tsin 2 # so that the condition 
for onset of decoherence (and loss of entanglement) is 
(277/7r 2 )( 7 2 /Aw)(P/V A )(T c Tsin 2 6l ) -1. Recent work 
\\?\ also found that noise will severely impact the perfor- 
mance of B((f>). 

We now show how a GQC containing noisy controlled- 
phase gates impacts the computational efficiency of the 
Shor algorithm for factoring an integer N |l4| . Num- 
ber theoretic arguments reduce factoring to finding the 
period r of the function F(a) — y a mod N, where y is 
co-prime with TV. The algorithm begins by preparing 
the GQC in the state |^ ) = (1/y/A+l) J2f =0 \jr + I)- 
Here r is the period of F(a); I is the result of a measure- 
ment carried out during preparation of IV'o); and A is the 
largest integer such that Ar + I < q, where q is chosen 
such that TV 2 < q < 2N 2 . The algorithm implements 
a discrete Fourier transform modulo q (DFT q ) on |^o) : 
l^i) = DFT q \^ ) = f(c)\c), where 



/(c) 



<1 



q/r-l 

E ex p 

3=0 



2m 



(jr + l)c + iT s (c) 



(8) 



The DFT q requires the application of L(L — l)/2 
controlled-phase gates, where L = log 2 q. Noisy con- 
trol of this gate causes the stochastic phase shift T s (c) to 
appear in eq. ©■ One can show 01 that T s (c) = 0, and 
the variance r 2 (c) is given by eq. (0 with i] = L(L — 1)/2; 

a 2 = P/VAuj; and I kj -> Y,Lo Ay The nnal ste P in the 
algorithm is a measurement which produces the result c 

with noise-averaged probability P{c) — |/(c)| 2 . In the 
absence of noise, constructive interference in eq. JSJ oc- 
curs when | re — c'q\ < r/2, which determines a unique 
c', and consequently, P(c') = P{c). The measurement 
result c yields the period r only if c' and r are co- 
prime. The success probability for the algorithm is then 
Psuc = ^Z' c i P(c'), where in the primed sum only those 
c' appear that are both less than, and co-prime with, r. 
One can show that: 



l/(c)| 2 



2r 



1 

q 

g/r-l 

EE- 

k=Q j>k 



Mi - k) 



<i 



(rc mod q) 



D(k,j), (9) 



where D(k,j) is the decoherence factor due to noise. As 

above, D(k,j) = exp — T 2 (k,j)/2 which, by eq. JSJ), 

goes smoothly to zero as the noise variance a 2 increases. 
In the decoherence limit only the first term on the RHS 

of eq. © survives so that |/(c)| 2 ~ l/q ~ 1/N 2 . This re- 
sult is independent of c so that P suc ~ (l/iV 2 )0(r), where 
4>{r) is Euler's Phi function which gives the number of in- 
tegers that are less than and co-prime with r. For large 



N, r < N and cj>(r) ~ r/logr - TV/ log TV. Thus, in the 
decoherence limit, P suc - l/(JVlogJV) = 2- lo e Ar /logA^. 
We see that the algorithm must be run on average 2 log N 
times to obtain P suc ~ 1. This is exponential in the 
problem size log N, indicating that sufficiently strong 
noise destroys the computational efficiency of the Shor 
algorithm. This contrasts with the (computationally ef- 
ficient) noiseless Shor algorithm whose runtime scales 
linearly with logiV. Onset of decoherence occurs when 
ff(c) - 4tt 2 , and for a GQC, £)f=o I l kj - t c T sin 2 O so 
that decoherence occurs when: 



VJ (aw) 



> 



TL(L - 1) 7 2 sin 2 



(10) 



Adiabatic operation of a GQC requires T to be large; and 
the desire to factor large N means L — log 2 q ~ log 2 N 
will also be large. Thus the RHS of eq. (|1U|) is expected 
to be small, so that managing noise will in fact be a 
significant issue for a GQC after all. Eq. flflj) allows a 
quantitative estimate of how much noise can be tolerated 
by a GQC before decoherence due to noisy control under- 
mines its operation. We believe that the above analysis is 
the first to quantitatively demonstrate how noisy control 
destroys the efficiency of Shor algorithm when run on a 
sufficiently noisy GQC. 
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